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REBREATHING METHODS INCLUDING OSCILLATING, 
SUBSTANTIALLY EQUAL REBREATHING AND NON-REBREATHING PERIODS 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[0001] The present invention relates generally to methods for noninvasively determining 
the pulmonary capillary blood flow ("PCBF") or cardiac output ("CO") of a patient. Particularly, 
the present invention relates to so-called differential Fick techniques for determining PCBF or CO 
including partial rebreathing techniques. 

Background of the Related Art 

[0002] Pulmonary capillary blood flow and cardiac output are examples of various 
hemodynamic parameters that may be monitored in critically ill patients. Cardiac output is the 
sum of blood flow through the lungs that participates in gas exchange, which is typically referred 
to as pulmonary capillary blood flow, and the blood flow that does not participate in gas 
exchange, which is typically referred to as intrapulmonary shunt flow or venous admixture. 

[0003] Conventionally, pulmonary capillary blood flow and cardiac output have been 
measured by direct, invasive techniques, such as by indicator dilution. Indicator dilution includes 
introducing a predetermined amount of an indicator into the bloodstream through the heart of a 
patient and analyzing blood downstream from the point of introduction to obtain a time vs. 
dilution curve. Thermodilution, in which room temperature or colder saline solution, which is 
also referred to as "cold" saline, is used as the indicator, is a widely employed type of indicator 
dilution. The cold saline is typically introduced into the right heart bloodstream of a patient 
through a Swan-Ganz catheter, which includes a thermistor at an end thereof The thermistor is 
employed to measure the temperature of the blood after it has passed through the right heart, or 
downstream from the point at which the cold saline is introduced. A thermodilution curve is then 
generated from the data, from which the cardiac output of the patient may be derived. 
Thermodilution and other indicator dilution techniques are, however, somewhat undesirable due 



to the potential for harm to the patient that is associated with inserting and maintaining such 
catheters in place. 

[0004] Less invasive indicator dilution methods that do not require that a catheter pass 
through the valves of the right side of the heart have also been developed. These less invasive 
methods include the so-called "transpulmonary indicator methods", which include the placement 
of probes in the esophagus or trachea (e.g., in Doppler/Transesophageal echocardiography). 
While the use of esophageal or endotracheal probes may seem less invasive than the introduction 
of a catheter into the heart of a patient, the potential for harming a patient exists nonetheless. 

[0005] Thus, safer, noninvasive techniques for determining pulmonary capillary blood 
flow and cardiac output have been developed. These noninvasive techniques are typically based 
on some form of the basic physiological principle known as the Fick principle: the rate of uptake 
of a substance by the blood or release of a substance from blood at the lung is equal to the blood 
flow past the lung and the content difference of the substance at each side of the lung. 

[0006] One variation of the Fick principle is the so-called carbon dioxide Fick equation: 
Qpcbf = Vco 2 /(CvC0 2 - CaC0 2 ), (1) 
where Q pcbf is pulmonary capillary blood flow, Vco 2 is carbon dioxide elimination, CvC0 2 is 
carbon dioxide content of the venous blood of the patient, and CaC0 2 is the carbon dioxide 
content of the arterial blood of the patient. 

[0007] Typically, a differential form of the carbon dioxide Fick equation is used to 
noninvasively determine the pulmonary capillary blood flow or cardiac output of a patient. Each 
of the differential Fick techniques for determining the pulmonary capillary blood flow or cardiac 
output of a patient are based on the fundamental premise that pulmonary capillary blood flow and 
cardiac output can be estimated based on the changes of other, measurable parameters when a 
change in the effective ventilation (i.e., the total ventilation less the wasted ventilation due to 
deadspace associated with the apparatus, the patient, or a combination thereof) occurs. When a 
differential form of the Fick equation is used, the pulmonary capillary blood flow or cardiac 



output of a patient may be determined on the basis of differences in Vco 2 , CaC0 2 , and CvC0 2 
between "normal" respiration and while a change in the effective ventilation of the patient is being 
induced. The following is an example of a differential Fick equation: 

QpcbfBD = YCQ 2 B - VCQ 2 D (2) 

(CvC0 2B - CvC0 2D ) - (CaC0 2B - CaC0 2D ), 
where Vco 2 B and Vco 2 D are, respectively, the carbon dioxide eliminations of the patient during 
"normal" breathing and while a change in the effective ventilation of the patient is being induced, 
CvC0 2B and CvC0 2D are the contents of carbon dioxide in the venous blood of the patient during 
the same periods, and CaC0 2B and CaC0 2D are the contents of carbon dioxide in the arterial 
blood of the patient during "normal" breathing and when the effective ventilation of the patient is 
changed, respectively. 

[0008] Typically, differential Fick techniques rely upon baseline measurements (i.e., 
taken during "normal" respiration) of Vco 2 and PetC0 2 . Once baseline data has been gathered, a 
change in the effective ventilation of the patient is induced. Once the Vco 2 and PetC0 2 values 
become stable with the change in effective ventilation, these parameters are again measured. The 
difference between the baseline values and those taken during the change in the effective 
ventilation of the patient are used to calculate the pulmonary capillary blood flow or cardiac 
output of the patient. When continually monitoring and updating the pulmonary capillary blood 
flow or cardiac output of a patient, a recovery period typically precedes reestablishment of the 
baseline values. The recovery period has been provided to facilitate the reestablishment of 
baseline levels prior to the start of rebreathing. Most investigators do not collect any data for 
analysis during the recovery period. 

[0009] The carbon dioxide Fick equation (1) and the differential Fick carbon dioxide 
equation (2) each require a determination of the Vco 2 of a patient. Carbon dioxide elimination is 
the net volume of carbon dioxide produced by the patient, or excreted from the body of a patient, 
during respiration. Therefore, carbon dioxide elimination is useful as an indicator of the metabolic 
rate of the patient. The Vco 2 of a patient may be noninvasively measured as the difference, per 
breath, between the volume of carbon dioxide inhaled during inspiration and the volume of carbon 



dioxide exhaled during expiration. Carbon dioxide elimination over a breath is typically calculated 
as follows: 

VC0 2 = / breath Vx/ C02 ^ (3) 

where V is the measured respiratory flow and f CQ2 is the substantially simultaneously detected 
carbon dioxide signal, or fraction of the respiratory gases that comprises carbon dioxide, or 
"carbon dioxide fraction". 

[0010] During rebreathing, the exhaled volume of carbon dioxide may change only 
slightly, while the inhaled volume of carbon dioxide, which is normally negligible, may increase 
substantially. As a consequence, the difference between the amounts of carbon dioxide that are 
exhaled and inhaled during rebreathing is reduced substantially, as is the carbon dioxide 
elimination of a patient. 

[0011] A determination of the CaC0 2 of a patient is typically based upon the measured 
y PetC0 2 of the patient. The PetC0 2 , after correcting for any deadspace in the patient's airway or 
CO in a ventilation circuit, is typically assumed to be approximately equal to the partial pressure of 

Lji. 

y carbon dioxide in the alveoli (PaC0 2 ) of the patient's lungs or, if there is no intrapulmonary 

I?. shunt, the partial pressure of carbon dioxide in the arterial blood of the patient (PaC0 2 ). Using a 

IH standard carbon dioxide dissociation curve, either the PetC0 2 measurement or the PaC0 2 

n calculation may be used to determine CaC0 2 . 

g : I 

l7\ [0012] When a change in the effective ventilation of a patient occurs, such as when a 

S3 patient inhales increased concentrations of carbon dioxide, CaC0 2 changes relatively quickly 
g spared to the rate of change in CvC0 2 , which has a higher carbon dioxide content CaCO, 
The content of carbon dioxide in the venous blood also changes relatively slowly because the 
body stores a large volume of carbon dioxide in other tissues. The carbon dioxide stores of an 
"average" human male may be as high as about 15 to 40 liters. Thus, the duration or magnitude 
of a change in the effective ventilation of a patient must be significant (e.g., the patient must inhale 
a significant amount of carbon dioxide) to effect a measurable change (e.g., increase) in the 
content of carbon dioxide in the patient's venous blood. Likewise, the Vco 2 of a patient changes 
at a faster rate than CvC0 2 . In fact, CvC0 2 changes so slowly relative to the rates at which 
CaC0 2 and the Vco 2 change that, when a change in effective ventilation occurs over a relatively 
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brief period of time (e.g., a few minutes or less), the CvC0 2 can be assumed to remain 
substantially the same (i.e., little or no change) over the time it takes to complete a conventional 
rebreathing maneuver. As the effects of rebreathing are delayed due to transport time and 
damping, conventional rebreathing processes typically employ a recovery period. 

[0013] Carbon dioxide elimination and the PetC0 2 are typically measured during both of 
the phases of a differential Fick technique. 

[0014] In one example of a known differential Fick technique for inducing a change in 
the effective ventilation of a patient, carbon dioxide may be added to the gases that are inhaled by 
the patient, either directly (e.g., by the addition of carbon dioxide from a cylinder or other external 
source) or by causing a patient to rebreathe previously exhaled gases. An exemplary differential 
Fick technique that has been employed, which is disclosed in Gedeon, A. et al. in 18 Med. & Biol 
Eng. & Comput. 411-418 (1980) (hereinafter "Gedeon"), employs a period of increased 
13 ventilation followed immediately by a period of decreased ventilation. When the technique 

disclosed in Gedeon or another so-called "rebreathing" process is used, the Vco 2 of the patient 

f ^ decreases to a level that is less than that which is measured during normal breathing. Rebreathing 

SLii 

ftj during which the Vco 2 decreases to near zero is typically referred to as "total rebreathing". 

It Rebreathing that causes some decrease, but not a total reduction of Vco 2 , is typically referred to 

;L as "partial rebreathing". These rebreathing processes may be used either to noninvasively 

ILy estimate the CvC0 2 , as in "total rebreathing", or to obviate the need to know CvC0 2 , as in 

*5 "partial rebreathing". 

-=* [0015] Rebreathing is typically conducted with a rebreathing circuit, which causes a 

Li 

patient to inhale a gas mixture that includes carbon dioxide. For example, the rebreathed air, 
which may be inhaled from a deadspace during rebreathing, includes air that has been exhaled by 
the patient (i.e., carbon dioxide-rich air). 

[0016] During total rebreathing, substantially all of the gas inhaled by the patient was 
expired during the previous breath. Thus, during total rebreathing, PetC0 2 is typically assumed to 
be equal or closely related to the partial pressure of carbon dioxide in the arterial (PaC0 2 ), venous 
(PvC0 2 ), and alveolar (PaC0 2 ) blood of the patient. Total rebreathing processes are based on the 
assumption that neither the pulmonary capillary blood flow or cardiac output, nor the CvC0 2 of 
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the patient, changes substantially during the rebreathing process. The partial pressure of carbon 
dioxide in blood may be converted to the content of carbon dioxide in blood by means of a carbon 
dioxide dissociation curve, where the change in the carbon dioxide content of the blood (CvC0 2 - 
CaC0 2 ) is equal to the slope(s) of the carbon dioxide dissociation curve multiplied by the 
measured change in PetC0 2 , as caused by a change in effective ventilation, such as rebreathing. 

[0017] In partial rebreathing, the patient inhales a mixture of "fresh" gases and gases 
that were exhaled during the previous breath. Thus, the patient does not inhale a volume of 
carbon dioxide as large as the volume of carbon dioxide that would be inhaled during a total 
rebreathing process. , 

[0018] As an example of a known partial rebreathing process, the NICO™ system 
offered by Novametrix Medical Systems Inc. of Wallingford, Connecticut, employs a 60 second 
baseline period, a 50 second rebreathing period, and/a 70 second recovery period. The complete 
rebreathing cycle lasts for about three minutes. Mother exemplary partial rebreathing process is 
disclosed in Capek, JM, and Roy, RJ, Noninvasive measurement of cardiac output using partial 
C0 2 rebreathing. IEEE Trans Biomed Eng l#88; 35:653-661. That rebreathing process has a 
total cycle time of about 3V4 minutes, withihe actual rebreathing phase lasting for about 30 
seconds. Gama de Abreu, M, et al., Partial carbon dioxide rebreathing: A reliable technique for 
noninvasive measurement of nonshunted pulmonary capillary blood flow, Crit. Care Med. 1997; 
25: 675-683, discloses a rebreathing'process with a 35 second rebreathing phase and a total cycle 
time, including baseline and recovery phases, of about 3 minutes. 

[0019] Conventional partial rebreathing processes typically employ a differential form of 
the carbon dioxide Fick equation, such as equation (2), to determine the pulmonary capillary 
blood flow or cardiac output of the patient without requiring knowledge of the carbon dioxide 
content of the venous blood of the patient since the carbon dioxide content of the venous blood of 
the patient is assumed to remain substantially the same (i.e., constant) in the periods during which 
measurements are obtained. 

[0020] Again, with a carbon dioxide dissociation curve, the measured partial pressure of 
end tidal carbon dioxide can be used to determine the change in content of carbon dioxide in the 
blood before and during the rebreathing process. Accordingly, the following equation can be used 



to determine pulmonary capillary blood flow or cardiac output when partial rebreathing is 
conducted: 

Q = AVC0 2 /sAPetC0 2 , (4) 
where s is the slope of the carbon dioxide dissociation curve. 

[0021] While partial rebreathing is the most commonly used method for causing a 
change in the effective ventilation of a patient, alternative differential Fick techniques for 
measuring pulmonary capillary blood flow or cardiac output have also been employed. Such 
differential Fick methods typically include a brief change of PetC0 2 and Vco 2 in response to a 
change in effective ventilation. This brief change can be accomplished by adjusting the respiratory 
rate, inspiratory and/or expiratory times, tidal volume, inspiratory pause, or positive-end 
expiratory pressure (PEEP) of the patient's respiration. 

[0022] While many existing differential Fick techniques provide reliable, noninvasively 
obtained measurements of pulmonary capillary blood flow or cardiac output, the lengths of time 
over which these techniques are effected are somewhat undesirable, especially in the common 
critical and intensive care situations in which it is desirable to more frequently update 
measurements of the pulmonary capillary blood flow or cardiac output of a patient. 

[0023] Accordingly, there is a need for a method of noninvasively calculating pulmonary 
capillary blood flow and cardiac output with increased frequency. 

BRIEF SUMMARY OF THE INVENTION 
[0024] The present invention includes a differential Fick technique for noninvasively 
determining the pulmonary capillary blood flow or cardiac output of a patient. The differential 
Fick method of the present invention includes two phases: a "normal" respiration phase and a 
phase in which a change in the effective ventilation of a patient is induced, which phase is referred 
to herein as a "change-inducing phase". These phases are abbreviated in duration relative to 
similar phases in known differential Fick techniques. The phases of the inventive differential Fick 
technique may be repeatedly cycled, or oscillated, with the durations of the normal respiration 
phase and the change-inducing phase being substantially the same. 
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[0025] In one aspect, the inventive differential Fick technique includes effecting a 
change-inducing phase in the respiration of a patient, allowing the respiration of the patient to 
return to normal, then immediately repeating the change-inducing phase of respiration. This 
method differs from conventional differential Fick techniques in that the typical recovery period, 
where a patient's respiration is allowed to return to normal, or baseline levels before again 
measuring respiratory carbon dioxide and flow is omitted. 

[0026] In another aspect of a differential Fick technique incorporating teachings of the 
present invention, the durations of the normal respiration and change-inducing phases are 
abbreviated relative to the time lengths of the corresponding phases in conventional differential 
Fick techniques. For example, each phase may have a duration of about 30 seconds. The length 
of an entire cycle of the differential Fick technique, measured as the difference in time between the 
end of one change-inducing phase and the end of another, immediately subsequent change- 
inducing phase is also shortened relative to the durations of conventional cycles of comparable 
differential Fick techniques. For example, a differential Fick technique conducted in accordance 
with teachings of the present invention may have a cycle time of about two minutes or less. 

[0027] Teachings of the present invention may be applied to rebreathing processes, as 
well as other differential Fick techniques. For example, changes in the respiration rate of a 
patient, changes in a patient's tidal volume, changes in the inspiratory pause of a patient, and 
changes in the positive-end expiratory pressure of a patient each induce changes in the effective 
ventilation of the patient that may be used in a differential form of the Fick equation to 
noninvasively determine the pulmonary capillary blood flow or cardiac output of the patient. 

[0028] Other features and advantages of the present invention will become apparent to 
those of ordinary skill in the art through consideration of the ensuing description, the 
accompanying drawings, and the appended claims. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 
[0029] In the drawings, which illustrate an exemplary embodiment for carrying out the 
invention: 
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[0030] FIG. 1 is a schematic representation which illustrates the componentry that may 
be utilized to measure respiratory profile parameters in a method incorporating teachings of the 
present invention; 

[0031] FIG. 2 A illustrates an idealized, bi-directional rebreathing cycle with Vco 2 values 
for different breaths depicted as diamonds and PetC0 2 values for various breaths shown as 
squares; 

[0032] FIG. 2B is a two-dimensional plot illustrating the use of a bi-directional 
rebreathing process to obtain three Vco 2 values and three values representative of the carbon 
dioxide content of the blood of a patient before, during, and after rebreathing; these three values 
have been used to substantially noninvasively determine the pulmonary capillary blood flow or 
cardiac output of the patient; 

[0033] FIG. 3 is a line graph that illustrates the Vco 2 of a patient during each of the 
before, during, and after phases of the bi-directional rebreathing method; 

[0034] FIG. 4 is a line graph that illustrates the PetC0 2 of a patient during each of the 
before, during, and after phases of the bi-directional rebreathing method; 

[0035] FIG. 5 A is a two-dimensional line graph illustrating a typical plot of Vco 2 on 
the y-axis and PetC0 2 on the x-axis; 

[0036] FIG. 5B is a two-dimensional line graph illustrating a plot of Vco 2 on the y-axis 
and CaC0 2 on the x-axis following modification of the Vco 2 and CaC0 2 data; 

[0037] FIG. 6 is a two-dimensional plot of a number of Vco 2 values against the same 
number of carbon dioxide content values obtained over a single bi-directional rebreathing cycle; 

[0038] FIG. 7 is an exemplary two-dimensional plot depicting modified Vco 2 and 
carbon dioxide content values from the same rebreathing cycle as that shown in FIG. 5 A; 

[0039] FIGs. 8A and 8B are two-dimensional plots illustrating an exemplary rebreathing 
technique and an accompanying method for modifying respiratory data to obtain an accurate best- 
fit line therethrough; 

[0040] FIG. 9A schematically depicts the timing of rebreathing and nonrebreathing 
phases of a conventional rebreathing process; and 
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[0041] FIG. 9B schematically depicts the timing of phases of an embodiment of a 
method incorporating teachings of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

[0042] As is known in the art, due to carbon dioxide buffering and large storage 
volumes of carbon dioxide in the body, the CvC0 2 of a patient changes very slowly relative to the 
rate at which the CaC0 2 of the patient varies when a change in the effective ventilation of the 
patient occurs. The inability of a signal, such as the CvC0 2 of a patient, to respond quickly to a 
change in effective ventilation is referred to as "damping" of the signal. 

[0043] FIG. 9A illustrates an example of a conventional rebreathing technique in 
which R indicates the periods during which rebreathing is effected and N designates the 
nonrebreathing periods, in which normal ventilation may occur. Conventional wisdom has 
indicated that differential Fick techniques should be short enough to avoid inducing changes in the 
CvC0 2 of a patient, but long enough to effect a substantial change in the CaC0 2 of the patient. 
Notwithstanding this conventional wisdom, it has been discovered by the inventors herein that any 
change in the effective ventilation of a patient sufficient to induce a change in the CaC0 2 of the 
patient without causing a significant increase in the CvC0 2 of the patient is useful. 

[0044] Accordingly, the present invention includes a method for inducing a change in 
the effective ventilation of a patient over a reduced period of time relative to that used in known 
rebreathing and other differential Fick techniques. An example of the method of the present 
invention is depicted in FIG. 9B. In FIG. 9B, R represents periods in which a change in the 
effective ventilation of a patient occurs, while N represents periods without a change in the 
effective ventilation of the patient, in which normal ventilation may occur. The change in effective 
ventilation occurs at an increased frequency relative to that at which changes in the effective 
ventilation of a patient are induced with previously known differential Fick techniques. 

[0045] Due to the decrease in total cycle time of the inventive differential Fick 
technique, as well as the increased portion of the cycle in which a change in effective ventilation is 
induced, an initial increase in the CvC0 2 of a patient may occur. Following such an initial 
increase in the measured CvC0 2 of the patient, however, further changes in the effective 
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ventilation of the patient, when conducted at a high enough frequency or oscillation rate, should 
not cause further, substantial changes in the CvC0 2 of the patient. Thus, in the differential Fick 
technique of the present invention, once the CvC0 2 of the patient reaches a new substantially 
steady state level, it remains at that level. Because CvC0 2 is damped, higher frequency changes, 
or oscillations, lead to a more stable CvC0 2 . 

[0046] In accordance with one aspect of the method of the present invention, a 
differential Fick technique may include a first phase, or ventilation state, in which a change in the 
effective ventilation is effected and a second phase, or ventilation state, in which "normal" 
ventilation, or breathing or respiration, occurs. 

[0047] Of course the change between the first and second ventilation states may be 
either abrupt, as depicted in FIG. 9B, or gradual. For example, when the change the effective 
ventilation of a patient is a change in the tidal volume of a patient, the change or oscillation 
between a minimum tidal volume and a maximum tidal volume may follow a sinusoidal curve 
(e.g., 400 ml, 410 ml, 420 ml, . . . 580 ml, 590 ml, 600 ml, 590 ml, 580 ml, . . . 420 ml, 410 ml, 
400 ml, 410 ml, 420 ml, . . .). 

[0048] The first and second phases, or ventilation states, of the inventive differential 
Fick technique may be effected for substantially the same amount of time, meaning that, while the 
first and second phases may have exactly the same duration, as depicted in FIG. 9B, one of the 
phases may alternatively be somewhat longer than the other. Either the first phase or the second 
phase of the differential Fick technique of the present invention may comprise as little as 
about 30% of the combined duration of the first and second phases or as much as about 70% of 
the combined duration of the first and second phases. For example, if the total cycle time of a 
differential Fick technique is about one minute, the first or second phase may last as little as 
about 18 seconds, while the other phase lasts about 42 seconds. Alternatively, the first and 
second phases may each last for about thirty seconds, as shown in FIG. 9B. 

[0049] The durations of each of the first and second phases, or ventilation states, may be 
optimized for a patient by evaluating the patient's ventilation (i.e., tidal volume x frequency). 
Optionally, the calculated pulmonary capillary blood flow or cardiac output of the patient may be 
evaluated along with the patient's ventilation to determine the optimal durations for each of the 
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first and second phases. For example, for higher cardiac output levels, it may desirable to reduce 
the time of each phase. 

[0050] As an alternative to oscillating between first and second phases, or ventilation 
states, respiratory data of a patient may be monitored before noninvasively measuring pulmonary 
capillary blood flow or cardiac output, prior to the initiation of a change in effective ventilation 
that facilitates such a noninvasive measurement (e.g., during a standby period). Accordingly, 
pulmonary capillary blood flow or cardiac output may be measured by initiating a change in the 
effective ventilation of a patient (i.e., the first phase) and using data that was obtained and stored 
during a standby period, prior to effecting a change in the effective ventilation of a patient. 

[0051] The cycle times of differential Fick techniques incorporating teachings of the 
present invention are preferably about two minutes or less. Cycle times of less than a minute are 
also within the scope of the present invention. 

[0052] Another aspect of the present invention is related to a discovery of the inventors 
herein that, in order to accurately, noninvasively measure pulmonary capillary blood flow or 
cardiac output, a recovery period is not necessary following the change in the effective ventilation 
of the patient. Accordingly, differential Fick techniques incorporating teachings of the present 
invention may lack the conventional recovery or stabilization period that typically follows 
inducement of a change in the effective ventilation of a patient. The first and second phases may 
be continuously cycled, one (e.g., the first phase) immediately following completion of the other , 
(e.g., the second phase). In addition, respiration (e.g., flow and carbon dioxide or oxygen levels) 
of the patient maybe continuously evaluated or monitored while the differential Fick technique of 
the present invention is being effected. Alternatively, one or more intermittent measurements may 
be obtained during each immediately sequential occurrence of the first and second phases. 

[0053] As an example of a differential Fick technique incorporating teachings of the 
present invention, partial rebreathing may be employed. In the partial rebreathing embodiment of 
the differential Fick technique, the first phase is a rebreathing phase, while the second phase is a 
non-rebreathing phase. In rebreathing, the Vco 2 of a patient is measured along with the PetC0 2 
of the patient or another indicator of the carbon dioxide content of the patient's blood (e.g., 
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CaC0 2 , pC0 2 a surrogate of PetC0 2 , such as the average pC0 2 over about the last 5% of the 
expired volume, etc.). 

[0054] FIG. 1 schematically illustrates an exemplary method of substantially 
noninvasively monitoring the respiration of a patient and of measuring the flow rates and carbon 
dioxide concentration of gas mixtures that are inhaled and exhaled by a patient 10 over the course 
of the patient's breathing, such as during normal respiration or during known rebreathing 
techniques. A respiratory circuit 1 1 is placed in communication with the airway of patient 10. A 
flow sensor 12 of a known type, such as the differential-pressure type respiratory flow sensors 
manufactured by Novametrix Medical Systems Inc. ("Novametrix") of Wallingford, Connecticut 
(e.g, the Pediatric/ Adult Flow Sensor (Catalog No. 6717) or the Neonatal Flow Sensor 
(Catalog No. 6718)), is positioned along ventilation circuit 11. Flow sensor 12 may be 
operatively attached to a ventilation apparatus (not shown) or patient breathing mask or 
mouthpiece, as well as respiratory flow sensors based on other operating principles and 
manufactured or marketed by others, may be used to measure the flow rates of the breathing of 
patient 10, 

[0055] A carbon dioxide sensor 14, such as the CAPNOSTAT® carbon dioxide sensor 
and a complementary airway adapter (e.g., the Pediatric/ Adult Single Patient Use Airway Adapter 
(Catalog No. 6063), the Pediatric/ Adult Reusable Airway Adapter (Catalog No. 7007), or the 
Neonatal/Pediatric Reusable Airway Adapter (Catalog No. 7053)), which are manufactured by 
Novametrix, as well as main stream and side stream carbon dioxide sensors manufactured or 
marketed by others, may be positioned along ventilation circuit 11 so as to measure the carbon 
dioxide concentration of gas mixtures that are inhaled and exhaled by patient 10. 

[0056] Flow sensor 12 and carbon dioxide sensor 14 are connected to a flow monitor 16 
and a carbon dioxide monitor 1 8, respectively, each of which may be operatively associated with a 
computer 20 so that data from the flow and carbon dioxide monitors 16 and 18 representative of 
the signals from each of flow sensor 12 and carbon dioxide sensor 14 may be detected by 
computer 20 and processed according to programming (e.g., by software) thereof. Preferably, 
raw flow and carbon dioxide signals from the flow monitor and carbon dioxide sensor are filtered 
to remove any significant artifacts. As respiratory flow and carbon dioxide pressure 
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measurements are made, the respiratory flow and carbon dioxide pressure data may be stored by 
computer 20. 

[0057] Each breath, or breathing cycle, of patient 10 may be delineated, as known in the 
art, such as by continuously monitoring the flow rate of the breathing of patient 10. 

[0058] In order to effect rebreathing, a deadspace 22, or carbon dioxide source, 
communicates with the airway of patient 10. During the non-rebreathing phase, communication 
between deadspace 22 and the airway of patient 10 is interrupted. 

[0059] In partial rebreathing in accordance with teachings of the present invention, a 
baseline may be established during the non-rebreathing phase, in which carbon dioxide elimination 
and the partial pressure of end tidal carbon dioxide are measured. The non-rebreathing phase is 
then immediately followed by a rebreathing phase, wherein a change in the CaC0 2 of the patient is 
induced and Vco 2 and PetC0 2 are again measured. The rebreathing phase is then immediately 
followed by another non-rebreathing phase, wherein Vco 2 and PetC0 2 are again measured. 

[0060] The differential Fick technique of the present invention may be used with 
conventional rebreathing maneuvers and processes, as well as other known rebreathing maneuvers 
and processes, which are modified by either shortening or completely eliminating the conventional 
recovery or stabilization periods of these maneuvers and processes. For example, the differential 
Fick technique of the present invention may be used with the so-called "bi-directional" process 
disclosed in U.S. patent application serial no. 09/150,136, filed September 9, 1998 (hereinafter 
"the '136 Application"), the disclosure of which is hereby incorporated in its entirety by this 
reference, or in the so-called "best-fit line" method, which is disclosed in U.S. patent application 
serial no. 09/510,702, filed on February 22, 2000 (hereinafter "the £ 702 Application"), the 
disclosure of which is hereby incorporated in its entirety by this reference. 

Bi-Directional Rebreathing 
[0061] In the bi-directional rebreathing process, as disclosed in the ' 136 Application, 
respiratory carbon dioxide and flow measurements are obtained in three phases: a "before" 
rebreathing phase, a "during" rebreathing phase, and an "after" rebreathing phase. When 
teachings of the present invention are applied to the bi-directional rebreathing method, 
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measurements obtained during a first non-rebreathing phase provide data for the "before" 
rebreathing period of a first rebreathing cycle, measurements obtained in the rebreathing phase 
provide data for the "during" rebreathing period of the first rebreathing cycle, and measurements 
obtained during the next non-rebreathing phase provide data for both the "after" rebreathing 
period of the first rebreathing cycle and the "before" rebreathing period of the next rebreathing 
cycle. 

[0062] In a first variation of the bi-directional rebreathing technique, which is useful 
when the CvC0 2 of a patient is changing, the rate at which the CvC0 2 changes is estimated, 
which may be useful for more accurately determining the pulmonary capillary blood flow or 
cardiac output of the patient. The rate of change in CvC0 2 may also be used to determine the 
rate or amount of change in the CvC0 2 of the patient between "normal" respiration and the 
change in effective ventilation. In a second variation of the bi-directional rebreathing technique, 
the pulmonary capillary blood flow or cardiac output of a patient may be determined without 
estimating the rate of change in CvC0 2 . The second variation of the bi-directional rebreathing 
technique is also useful for determining the pulmonary capillary blood flow or cardiac output of a 
patient when either CvC0 2 or cardiac output changes during the re-breathing process. 

[0063] Derivation of Formulae Employed in the Bi-directional Rebreathing Technique 
[0064] A differential form of the carbon dioxide Fick equation, similar to that employed 

in conventional partial re-breathing techniques, which is based on Vco 2 and CvC0 2 during 

"normal" respiration (B) and during the rebreathing process (Q), follows: 

Qbd = YC0 2 B - VCQ 2 D 

(CvC0 2 B - CvC0 2 D ) - (CaC0 2 b - CaC0 2 d ) (5) 

or 

Qbd = AVC0 2 BD 

ACvC0 2 BD - ACAC0 2 bd . (6) 
[0065] Similarly, another differential form of the carbon dioxide Fick equation, which is 

based on carbon dioxide elimination and carbon dioxide content measurements made during the 

rebreathing process and after rebreathing, which may also be employed to determine the 

pulmonary capillary blood flow or cardiac output of a patient, follows: 



16 



Qda = Yco 2 p - Vco 



■2 A- 



(CvC0 2 D - CvC0 2 J - (CaC0 2 d - CaC0 2 J (7) 



or 



Qda = AVCQ 



2 DA- 



ACvC0 2DA - ACaC0 2DA , (8) 

where CvC0 2 A is the CvC0 2 of the patient after rebreathing (A), or in the "after" phase. 

[0066] The two preceding differential forms of the carbon dioxide Fick equation may be 
combined to yield the following differential form of the carbon dioxide Fick equation: 
Q = AVCQ 2BD - AVco 2DA 



(ACvC0 2 BD - ACvC0 2 DA ) - (ACAC0 2 bd - ACaC0 2 da ). (9) 

[0067] Since CvC0 2 may change over tipfe, an accurate non-invasive Fick-based 
^^determination of the pulmonary capillary blocui'rlow or cardiac output of a patient may include an 
7J//estimation of the rate at which CvC0 2 ch^rlges. With an exemplary assumption that changes in 
^ / CvC0 2 are substantially linear over tl^rebreathing cycle and, therefore, that the rate of change is 
/ constant, the rate of change in CyC0 2 , represented as "k", may be determined with the following 



equation: 



k = ACvgq (10) 

At. 



Alternatively, the change in carbon dioxide content of the venous blood may be assumed to 
substantially follow a curoe of some other shape that is reasonably based on the character of the 
change in carbon dioxide content, such as an exponential curve, wherein the rate of change would 
also be exponential, or the curvesof a polynomial. As another alternative, the rate of change in 
CvC0 2 may be approximated by an artificial neural network or a radial basis function, as known in 
the art. 

[0068] When the change in CvCO^s assumed to be linear with respect to time and, 
therefore, the rate of change of CvC0 2 is assumed to be constant, the change in CvC0 2 between 
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the "before" and "during" phases and between the "during" and "after" phases can be expressed 
by the following equations: 

ACvC0 2BD = k(t B -t D ) (11) 

and 

ACvCO^^kCto-U (12) 
where t D , t B and t A represent the times at which the "before", "during" and "after" phases 
respectively occur. 

[0069] The foregoing equations for the change in CvC0 2 may be substituted into the 
differential form of the carbon dioxide Fick equation that considers the breathing of a patient 
during each of the "before", "during" and "after" phases and the "A" terms expanded to yield the 
following form of the carbon dioxide Fick equation, which accounts for any changes in CvC0 2 
and is, therefore, useful in the bi-directional rebreathing method: 

Q = VCQ 2 B _±VC0 2 A - 2 ' VCQ 2 D 

k • (t B + t A - 2 * t D ) - (CaC0 2B + CaC0 2 a - 2 • CaC0 2D ). (13) 

If, however, t D -t B =t A -t D as is probable in the differential Fick technique of the present invention, 
then t A +t B =2-t D , and it would not be necessary to calculate k, as k would be multiplied by zero. 
Accordingly, if t D -t B =t A -t D , such as when the durations of the first and second phases are the same, 
the following equation could be employed to determine the pulmonary capillary blood flow of a 
patient: 

Q =__Vco 2 P + Vco 2 A - 2 ■ Vco 2 D 

- (CaC0 2 b + CaC0 2 A - 2 • CaC0 2 D ). (14) 

[0070] Due to the assumption that the pulmonary capillary blood flow and cardiac 
output of a patient remain substantially constant from the "before" phase to the "after" phase, the 
differential carbon dioxide Fick equations for determining pulmonary capillary blood flow or 
cardiac output over the "before" and "during" phases (Q BD ) and for determining cardiac output 
over the "during" and "after" phases (Q DA ) may be. employed to estimate k, the rate of change in 
CvC0 2 , as follows: 

Qbd = Qda, (15) 

thus, 
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AVC0 2BD = AVco 2D> , (16) 



ACvC0 2 BD - ACaC0 2 bd, ACvC0 2 da - ACaC0 2 da , 

which may be rearranged as: 

A Vco 2 BD • ACvC0 2 DA - AVco 2 DA * ACvC0 2 BD - 

AVco 2BD • ACAC0 2DA - AVco 2DA • ACaC0 2BD . (17) 

The equations for ACvC0 2BD (11) and ACvC0 2DA (12) are then substituted into the preceding 
equation (17) to yield the following equation: 

AVco 2BD • k(t D - tj - AVco 2DA • k(t B - t D ) - 

AVco 2 BD • ACaC0 2 da - AVC0 2 DA • ACaC0 2 bd , (1 8) 

which may be rearranged to provide the following equation for k, the rate of change in Vco 2 : 

k = AVco 2 BD - ACaCQ 2 d a - AVco 2 DA ■ ACaCQ 2 bd 

AVC0 2BD (t D - tj - AVC0 2DA (t B - t D ). (19) 

[0071] Use of the Bi-Directional Rebreathing Technique While Cardiac Output is 
Changing to Noninvasively Determine Pulmonary Capillary Blood 

[0072] Equation (14) above may also used to determine the pulmonary capillary blood 
flow or cardiac output of a patient if the pulmonary capillary blood flow or cardiac output of the 
patient changes during the bi-directional rebreathing process. This can be shown by assuming that 
CvC0 2 does not change during re-breathing: 

CvC0 2 = CvC0 2 B = CvC0 2 D = CvC0 2 A . (20) 

[0073] The Fick equation can be used to express Vco 2 in terms of CvC0 2 and CaC0 2 or 
CaC0 2 during each of the before, during, and after phases of the bi-directional partial rebreathing 
method: 

Vco 2B = Q B (CvC0 2B - CaC0 2B ); (21) 
Vco 2 D = Q D (CvC0 2 D - CaC0 2 d ); and (22) 
Vco 2 A = Q A (CvC0 2 A - CaC0 2 J. (23) 
Substituting equations (21)-(23) for the Vco 2B , Vco 2D , and Vco 2A terms of equation (14) along 
with the assumption that CvC0 2 is not changing during rebreathing in equation (20) provides the 
following equation: 

Q = Q P (CvC0 2 - CaCCX p ) + QpfCvCO, - CaCCX p )- 2Q a (CvCQ 2 - CaCO ? ^ 
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- (CaC0 2B + CaC0 2A - 2 • CaC0 2D ). 



(24) 



[0074] When it is assumed that CvC0 2 does not change during rebreathing, it may be 

implied that once the rebreathing process is stopped, CaC0 2 and CaC0 2 return to substantially 

the same levels of these parameters prior to the rebreathing process, or the same levels of these 

parameters if rebreathing had not been conducted or another change in the ventilation of a patient 

had not been induced: 

CaC0 2 b = CAC0 2 a = CaC0 2 ^ (25) 

where CaC0 2 ^ is the content of carbon dioxide in the capillaries surrounding the alveoli during 

time periods in which a change in the ventilation of a patient has not been induced (e.g., as is 

induced during rebreathing), such as during the second phase of the differential Fick technique of 

the present invention. Thus, CaC0 2NR may be substituted for both CaC0 2B and CaC0 2 A in 

equation (24), providing the following equation: 

Q = O p CvCO ,-O p 'CaC0 2 >T p+O a 'CvCO ,-O a CaCO o >T p-20 r CvC0 2 ±Q D CaCO . D 

-2(CaC0 2 - CaC0 2 d ). (26) 

By factoring certain terms, equation (26) can be rewritten as follows: 

Q = rO P + O, -20 ,0 - CvCXX - (O, + 0 P ^ • CACO , >T + 20 p -CaCO, p 

-2(CaC0 2 ^ - CaC0 2 d ). (27) 

[0075] If it is assumed that the pulmonary capillary blood flow or the cardiac output of 

the patient changes linearly over time, or that the difference in pulmonary capillary blood flow or 

cardiac output between before and during phases of the bi-directional partial rebreathing method 

is equal to the difference in pulmonary capillary blood flow or cardiac output between the during 

and after phases, then the pulmonary capillary blood flow or cardiac output of the patient in the 

during phase of rebreathing may be expressed by the following equation: 

Qd - '/ 2 (Q A + Q B ). (28) 

Substituting equation (28) into equation (27) provides the following equation: 

Q = [Qa + Qb - (Qa + Op^- CvCO, - (0, + O pVCACO , h> + 2CACO , P (OA + QB) 

-2(CaC0 2NR -CaC0 2D ), (29) 

which is the equivalent of the following equation: 
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Q = V4 (Q A + Qb), 



(30) 



which is equal to Q D . 

[0076] Thus, when the pulmonary capillary blood flow or cardiac output of a patient 
changes at a constant rate during the bi-directional partial rebreathing method, the cardiac output 
or pulmonary capillary blood flow of a patient is equal to the pulmonary capillary blood flow or 
cardiac output of the patient as measured in the during phase of partial rebreathing. Accordingly, 
equations (13) and (14) may be used with the bi-directional partial rebreathing process to 
accurately and noninvasively determine the pulmonary capillary blood flow or the cardiac output 
of a patient while the pulmonary capillary blood flow or cardiac output of the patient is changing. 

[0077] The bi-directional partial rebreathing method, as embodied in equations (13) and 
(14), may also be useful for noninvasively determining the pulmonary capillary blood flow or 
cardiac output of a patient while both the CvC0 2 and the pulmonary capillary blood flow or 
cardiac output of the patient are changing. 

[0078] In practicing the bi-directional rebreathing method, a system such as that 
described in reference to FIG. 1 is used and the patient's breathing is monitored during 
rebreathing and non-rebreathing phases to detect the amount of C0 2 exhaled by the patient and 
the flow rate of the patient's respiration during these phases, from which Vco 2 and CvC0 2 may be 
determined. 

[0079] The graph of FIG. 2 A illustrates the various measurements that may be obtained 
during bi-directional rebreathing. As illustrated, the typical changes in the VC0 2 (shown as 
diamonds) and carbon dioxide content measurements (e.g., PetC0 2 , shown as squares) that may 
occur between the baseline breathing (i.e., "before" rebreathing), "during" rebreathing, and the 
recovery or stabilization (i.e., "after" rebreathing) periods of an idealized (i.e., without noise) 
bi-directional rebreathing cycle. During rebreathing, Vco 2 changes from a baseline value (e.g., 
about 200 ml/min) to a during rebreathing plateau (e.g., of about 100 ml/min.) within about 3 or 4 
breaths, whereas carbon dioxide content may take longer to change from a baseline value 
(e.g., 35 mmHg) to a plateau (e.g., about 39 mmHg). 

[0080] FIG. 2B is a two-dimensional plot that illustrates that single values, the plateau 
values, from each of the before, during, and after rebreathing phases of a bi-directional 
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rebreathing process, such as that illustrated in FIG. 2A, are used to estimate pulmonary capillary 
blood flow or cardiac output. 

[0081] The difference between the volume of carbon dioxide exhaled and the volume of 
carbon dioxide inhaled by a patient, which may be used to estimate the Vco 2 of the patient, is 
determined before (Vco 2 B ), during (Vco 2D ), and after (VC0 2 J rebreathing. FIG. 3 is a graph 
that illustrates Vco 2 during each of the before, during, and after phases of the rebreathing process 
of the present invention. 

[0082] The PetC0 2 of the patient is also m4sured for each of the "before", "during" 
and "after" phases. As PetC0 2 , when corrected for parallel deadspace (of non-perfused alveoli), 
is assumed to be equal to the partial pressure of carbon dioxide in the alveolar blood (PaC0 2 ) and 
the partial pressure of C0 2 in the arterial blood (PaC0 2 ), a carbon dioxide dissociation curve may 
be employed with the end tidal carbon dioxide partial pressure measurements, as known in the art, 
to determine the content of carbon dioxdde in blood of the alveoli (CaC0 2 ) of the lungs of the 
patient that participate in the exchange of blood gases, which alveoli are typically referred to as 
"perfused" alveoli, for each of the/before, during, and after rebreathing phases. CaC0 2 is 
assumed to be equal to the content of carbon dioxide in the arterial blood (CaC0 2 ). FIG. 4 is a 
graph that illustrates the etCG) 2 measured during each of the before, during, and after phases of 
the rebreathing process of me present invention. 

[0083] Determining Pulmonary Capillary Blood Flow or Cardiac Output 
[0084] In determining the pulmonary capillary blood flow or cardiac output of a patient 
when CvC0 2 changes, the differences between Vco 2 before rebreathing and during rebreathing, 
which difference is also referred to as "AVco 2BD ", and during rebreathing and after rebreathing, 
which difference is also referred to as "AVC0 2DA ", are determined. The differences between the 
CaC0 2 before and during rebreathing, which difference is also referred to as "ACaC0 2BD ", and 
during rebreathing and after rebreathing, which difference is also referred to as "ACaC0 2DA ", are 
also determined. 

[0085] These differences may then be used to calculate the rate at which the content of 
carbon dioxide in the venous blood of the patient changes. An exemplary equation for estimating 
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the rate of change in the content of carbon dioxide in the patient's venous blood (k), which 
assumes that the change is linear with time and, therefore, that the rate of change is constant, 
follows: 

k = AVco 2 PP • ACaCO, da - AVco 2 DA - ACaCQ 2 bd (31) 
AVco 2BD (t D - tj - AVco 2 

DA i^B " tb). 

[0086] Once the rate of change in the content of carbon dioxide in the patient's blood 

has been estimated, the pulmonary capillary blood flow or cardiac output of the patient may be 

accurately determined, as follows: 

Q = VC0 2 b_±VC0 2 A - 2VCQ 2 D (32) 

k ■ (t B + t A - 2t D ) - (CaC0 2 b + CaC0 2 a - 2CaC0 2 d ) 

[0087] Alternatively, the times and constant may be omitted from the previous equation 

and pulmonary capillary blood flow or cardiac output of the patient determined by use of the 

following equation, which is useful when either CvC0 2 or pulmonary capillary blood flow or 

cardiac output changes during re-breathing: 

Q = __Vco 2 B + VCQ 2 A - 2VCQ 2 D (33) 

- (CaC0 2B + CaC0 2a - 2CaC0 2D ). 



Best-Fit Line Method 

[0088] By way of contrast with the use of measurements at the plateaus of each of the 
phases, as depicted in FIGs. 2A and 2B, in conventional rebreathing processes and the 
bi-directional rebreathing process, Vco 2 and carbon dioxide content data are continually measured 
in the rebreathing method of the '702 Application. As a result, a plot of the measurements may 
have the appearance of the graph shown in FIG. 5 A, with data at 100 being based on before 
rebreathing measurements, data along arrow 102 being based on during rebreathing 
measurements, and data along arrow 104 being based on after rebreathing measurements. These 
data may be obtained by use of a single rebreathing cycle, over the course of a number of 
rebreathing cycles, at one or more discrete time intervals, or on a breath-by-breath basis, where 
data is continually measured, calculated, and analyzed in accordance with the method of the 
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invention so as to continually update or monitor the pulmonary capillary blood flow or cardiac 
output of a patient. 

[0089] The best-fit line method also includes use of a differential form of the carbon 

dioxide Fick equation to calculate pulmonary capillary blood flow or cardiac output as the ratio of 

a change in carbon dioxide elimination, or Vco 2 , to a change in the content of carbon dioxide, or 

CaC0 2 , in the arterial blood of a patient: 

Q = AVCO, (34) 
ACaC0 2 . 

As explained previously herein, CaC0 2 can be noninvasively estimated by determining the 
PetC0 2 . PetC0 2 may be converted to CaC0 2 by use of a standard carbon dioxide dissociation 
curve, as is known in the art, by use of the following equation: 

ACaC0 2 = sAPetC0 2 , (35) 
where s is the slope of the carbon dioxide dissociation curve and APetC0 2 is a change in the end 
tidal partial pressure of carbon dioxide of a patient effected by a change in ventilation. Thus, 
pulmonary capillary blood flow or cardiac output can also be calculated as follows: 

Q = AVco 2 /sAPetC0 2 . (36) 

[0090] Other indicators of the carbon dioxide content in the blood of a patient, such as 
pC0 2 , may be used in place of PetC0 2 or CaC0 2 to determine the pulmonary capillary blood flow 
or cardiac output of a patient. 

[0091] Once respiratory carbon dioxide pressure and flow measurements have been 
made, as depicted in FIG. 1, during both the first (e.g., rebreathing) and second (e.g., non- 
rebreathing) phases, these respiratory carbon dioxide pressure and flow data are used, as known 
in the art, to calculate Vco 2 and PetC0 2 , as well as the changes in Vco 2 and PetC0 2 that occur 
with the change in effective ventilation. 

[0092] The calculated Vco 2 and PetC0 2 data are then used to determine the pulmonary 
capillary blood flow or cardiac output of the patient, such as by use of any of the various Fick 
equations presented above. 

[0093] As an alternative, the pulmonary capillary blood flow or cardiac output of a 
patient can be determined over the course of a plurality of breaths by expressing the calculated 
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Vco 2 data and CaC0 2 data or data of another indicator of the content of carbon dioxide in the 
blood of a patient, such as PetC0 2 or pC0 2 , in two dimensions, such as on a two-dimensional 
(X-Y) line graph, with Vco 2 data points being measured on the y-axis and PetC0 2 data points 
being measured on the x-axis, then identifying a line that best fits the data, which is also referred 
to herein as a best-fit line. 

[0094] For example, the equation for the best fit line is: 

y = mx + b (37) 

or 

m = y-b (38) 

where y is the y-axis ordinate of a data point, x is the x-axis ordinate of the same data point, m is 

the slope of the line, and b is the offset value for the line. If Vco 2 is measured on the y-axis and 

CaC0 2 is measured on the x-axis, then 

m = Vco 2 - b (39) 
CaC0 2 . 

The negative slope (i.e., -1 x m) of the best-fit line through the Vco 2 -CaC0 2 data would be equal 
to the pulmonary capillary blood flow or cardiac output of the patient: 

-m = Q. (40) 
[0095] The best-fit line for the VC0 2 and CaC0 2 data is preferably determined by use of 
known linear regression techniques or any other known methodology for determining the 
relationship between two variables. The method of linear regression provides an accurate 
pulmonary capillary blood flow or cardiac output value based on a large number of Vco 2 and 
CaC0 2 data obtained over the course of one or more changes in effective ventilation. When linear 
regression is used, the slope (m) of the best-fit line for the data is calculated as follows: 

m = Lxy/Lxx (41) 
and the offset (b) of the line is calculated by the following equation: 

b - Ey/n - m x Ex/n, (42) 

where 

Lxx = Ex 2 - (Ex x ix)/n, (43) 
Lyy = Ey 2 - (Ey x Ey)/n, and (44) 
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Lxy = Exy - (Ex x Ey)/n, (45) 
and where n is the number of data points in the plot, Ex is the sum of all x-ordinate (i.e., CaC0 2 
content) values, Ey is the sum of all y-ordinate (i.e., VC0 2 ) values, Ex 2 is the sum of the square of 
all x-ordinate values, Ey 2 is the sum of the square of all y-ordinate values, and Sxy is the sum all 
paired x- and y-ordinate values multiplied by each other. 

[0096] When linear regression is used to determine the location and orientation of a 
best-fit line, a correlation coefficient (r) that quantifies the accuracy with which the best-fit line 
correlates to the Vco 2 and CaC0 2 data can also be calculated as follows: 

r = (Lxy x Lxy)/(Lyy x Lxx). (46) 
Alternatively, any other measure of the quality of fit that quantifies the accuracy with which the 
best-fit line correlates to the Vco 2 and CaC0 2 data may be used. 

[0097] Correlation coefficients range from 0 to 1.0, where a correlation coefficient of 0 
indicates that no linear correlation exists between the x-ordinate and the y-ordinate data and a 
correlation coefficient of 1.0 indicates that the x-ordinate and y-ordinate data are perfectly linearly 
correlated (i.e., all of the Vco 2 -CaC0 2 data points are located on the same straight line). 

[0098] The Vco 2 -CaC0 2 data points measured before and during rebreathing, however, 
are rarely located on the same straight line. One reason for this is that, during rebreathing 
maneuvers, the Vco 2 signal typically leads the PetC0 2 signal and, thus, the CaC0 2 by about one 
breath. In addition, Vco 2 is calculated on the basis of signal components that have higher 
frequencies than do the PetC0 2 signal. As a result, when the Vco 2 and CaC0 2 measurements 
calculated over a period of time are plotted against one another on a two-dimensional (X- Y) line 
graph, the result typically appears as an arc or a loop, as shown in FIGs. 5 A and 6, rather than as 
a straight line, depending on the amount of data calculated and the duration of rebreathing. 
Moreover, Vco 2 and CaC0 2 measurements may be calculated on the basis of respiratory flow and 
carbon dioxide pressure data obtained during spurious breaths. Such data do not relate to the 
pulmonary capillary blood flow or cardiac output measurement. Vco 2 and CaC0 2 calculations 
that are based upon such spurious data act as noise that may result in miscalculation of a best-fit 
line through the calculated VC0 2 and CaC0 2 data. As a result, the correlation coefficient of a 
best-fit line to the data is typically much less than 1.0. 
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[0099] The measured respiratory flow and carbon dioxide pressure data or the 
calculated Vco 2 and CaC0 2 data can be modified to increase the correlation coefficient between 
the Vco 2 and CaC0 2 data and the best-fit line therefor. Preferably, a linear transform is used to 
increase the correlation coefficient. A linear transform may be used to delay the calculation of a 
Vco 2 data point to accurately coincide therewith a CaC0 2 data point based on measurements 
taken during the same breath. The measured or calculated data may also be filtered by use of a 
linear transform. 

[0100] In an exemplary method for increasing the correlation coefficient between the 
Vco 2 and CaC0 2 data and the best-fit line therefor, a filter is applied to the calculated Vco 2 or 
CaC0 2 data. Known analog or digital low-pass, high-pass, or band pass filters, including adaptive 
filters, may be employed. Linear or nonlinear filters may be used. Preferably, a first order (single 
pole) infinite impulse response (IIR) digital filter is employed to filter the VC0 2 calculations in a 
manner that improves the correlation between the Vco 2 calculation and the lagging 
PetC0 2 /CaC0 2 calculation. The equation for such a filter is: 

Vco 2 >] - a x VC(V[n-l] + Q-a) * VC0 2 [n], (47) 
where Vco 2 [n] is the most recently calculated, unfiltered Vco 2 data point, Vco 2 '[n-l] is the 
previous, filtered Vco 2 data point, Vco 2 "[n] is the new "filtered" value based on Vco 2 [n] and 
obtained by use of the filter, and a is the filter coefficient. The filter coefficient, a, has a range 
of 0 to 1.0. The greater the value of a, the more profoundly the most recently calculated data 
point is filtered and, conversely, the lower a values cause the most recently calculated data points 
to be filtered to a lesser degree. When a is equal to zero, the most recently calculated data point 
is not filtered. 

[0101] Due to anatomical and physiological differences between different patients, 
different patients have differing optimal filter coefficients, a. In addition, as anatomical and 
physiological changes may occur in a patient over time, the optimum filter coefficient, a, to be 
used in filtering the Vco 2 or CaC0 2 values calculated from the patient's breathing may also vary 
over time. Any known optimization method or search algorithm may be employed to select 
optimal filter coefficient, a. 
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[0102] As an example of one way in which an optimal filter coefficient may be selected, 
a is first set to a default value (e.g., 0.85) and the calculated Vco 2 or CaC0 2 values are filtered on 
the basis of the default filter coefficient, a. The linear regression is then performed to obtain a 
best-fit line. If the correlation coefficient of the best-fit line calculated with the just-filtered data is 
less than the correlation coefficient of the immediately preceding best-fit line, which was 
calculated with unfiltered data or with a prior filter coefficient, then a predetermined a adjustment 
value (e.g., 0.01) is changed by multiplying the a adjustment value by -1 and by modifying the 
filter coefficient by adding the modified a adjustment value thereto. Otherwise, the filter 
coefficient, a, is modified by adding the unmodified a adjustment value thereto. The process of 
filtering the data based on a modified filter coefficient, obtaining a best-fit line for the data, 
comparing the correlation coefficient of the best-fit line to the correlation coefficient of the 
previous best-fit line, and adjusting the filter coefficient accordingly is then repeated a 
predetermined number of times (e.g., 50 times). The best-fit line with the greatest correlation 
coefficient, based on the unfiltered data and each set of filtered data, is selected to calculate the 
pulmonary capillary blood flow or cardiac output of the patient. When filtering is used, the Vco 2 - 
CaC0 2 plot preferably narrows, as depicted in FIGs. 5B and 7, to thereby increase the accuracy 
with which the location and orientation of a best-fit line can be established and, thus, to increase 
the accuracy of a pulmonary capillary blood flow or cardiac output determination based on the 
data. 

[0103] Another example of a method for increasing the correlation coefficient between 
the Vco 2 and CaC0 2 data and the best-fit line therefor, which is referred to herein as "clustering", 
includes the selection of data points that are grouped closely together. That is, the data points 
that are selected include those data points having a number of other data points within a 
predetermined range thereof. Data points that are not clustered are probably inaccurate or based 
on measurements taken during spurious breaths. As an accurate best-fit line through the data 
would likely be based on the clustered data, the data points that are not located in a cluster are not 
used in calculating the location and orientation of a best fit line for the data. 

[0104] Clustering of the data points may include normalization or transformation of the 
data such that ranges of the x-ordinate data (e.g., the CaC0 2 data) and the y-ordinate data (e.g., 
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the Vco 2 data) are substantially the same. Without such normalization, the data group (e.g., the 
Vco 2 data or the CaC0 2 data) with the highest range would dominate; the other data group 
would be less significant. 

[0105] An exemplary manner in which the data may be normalized includes use of the 
following normalization: 

x = (x-x)/o x , (48) 

where: 

x is the raw value, x is the mean value of all x-axis (e.g., CaC0 2 ) data in the plot, and a x is the 
standard deviation of all x-axis data in the plot. This normalization is applied to all x-axis values. 
A similar normalization scheme is applied to all of the y-axis values. 

[0106] The normalized data may then be clustered by searching for a predetermined 
number (e.g., 5) of the closest data points (e.g., Vco 2 or CaC0 2 data points) to each of the data 
points in a group. The differences between the analyzed data point and each of the predetermined 
number of closest data points are then added together and compared to a predetermined 
threshold. If the sum of the differences exceeds the predetermined threshold, the analyzed data 
point is discarded. Of course, the use of other clustering techniques to identify the most accurate 
data and to disregard probable inaccurate data are also within the scope of the present invention. 

[0107] Once clustering has been performed, the inverse of the normalization is 
calculated, or the normalization is undone, to provide an accurate determination of pulmonary 
capillary blood flow or cardiac output. An example of the manner in which the inverse of the 
normalization may be calculated includes use of the following equation: 

x = xo x + x. (49) 
This inverse of the normalization is applied to all of the clustered x-axis (e.g., CaC0 2 ) values. A 
similar inverse normalization scheme is applied to all of the clustered y-axis data. 

[0108] Clustering is one of many known techniques for determining outliers. Other 
known techniques for determining outliers may also be used in the method of the present 
invention. 

[0109] Alternatively, or in addition to disregarding probable inaccurate data points, in 
order to enhance the accuracy of the data, clustering can be used add synthetic data points. 
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Synthetic data points may be added to increase the cordon coefficient of the best-fit line to the 
data points on which the best-fit line is based. / 

[0110] Another exemplary method^fmodifying data in the best fit line method is 
depicted in FIGs. 8A and 8B. As witMhe^ltering and clustering methods for modifying data, the 

ethod depicted in FIGs. 8A ancLJJBmcludes selection of data points that are most likely to 
facilitate an accurate determj^n of the location and orientation of a best-fit line and, thus, of 
the pulmonary capillary^d flow or cardiac output of a patient. This method for modifying 
data includes itera^y examining data points and the distribution of the remaining data points 
relative to theytwo lines representing the range of possible pulmonary capillary blood flow 
measurements. 

[0111] As shown in FIGs. 8A and 8B, a line or the equation for a line 110 representing a 
minimum expected pulmonary capillary blood flow (i.e., -m^ PCBF..J and a line or the 
equation for a line 120 representing a maximum expected pulmonary capillary blood flow (i.e., - 
mune 12 o = PCBF^J are positioned to intersect at a data point 130. For example, when the x- 
ordirlte is based on CaC0 2 , line 110 may have a slope of -0.5, which represents a minimum 
expected pulmonary capillary blood flow of 0.5 L/min, and line 120 may have a slope of -20, 
which represents a maximum pulmonary capillary blood flow of 20 L/min. Of course, other 
pulmonary capillary blood flow values for lines 1 10 and 120 may also be used. 

[0112] Next, the number of other data points 130 located between lines 1 10 and 120 is 
determined. If the number of data points 130 beWen lines 110 and 120 is equal to or exceeds a 
threshold number, the analyzed data point 130 is/etained for a subsequent determination of the 
location and orientation of a best-fit line throu/h the data. Otherwise, the analyzed data point 130 
is discarded. The threshold number of data/oints that must be located between line 1 10 and 
line 120 for an analyzed data point to be/tained may be a predetermined value or determined by 
other means. As example, the threshol/number may be set to the median number of data points 
that are located between line 1 10 andCe 120 when each data point 130 of a set of data 
points 130 has been evaluated in accordance with the present embodiment of the method for 
modifying data. This process is Repeated until each data point 130 in a set of data points 130 has 
been so evaluated. FIG. 8A d/icts use of the data modification method on a data point 130 that 
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ill be retained, while FIG. 8B illustrates use of the present embodiment of the data modification 
method on another data pointp30' that will not be retained. 

[0113] FIGs. 5 A and 5B and FIGs. 6 and 7 illustrate the effect of modifying data to 
increase the accuracy with which the location and orientation of a best-fit line through the data 
may be determined. FIG. 6 illustrates a typical Vco 2 vs. CaC0 2 plot without such modification, 
where the plot appears as a loop. By way of contrast, FIG. 7 illustrates the closeness of the data 
when one or more of the embodiments of the method of the present invention are used to modify 
the data. FIGs. 5 A and 5B illustrate plots of Vco 2 and PetC0 2 data before and after modification, 
respectively. The increased closeness of the data points makes it possible to determine the 
orientation and location of a best-fit line therethrough with increased accuracy. 

[0114] Once all of the data points have been examined, the location and orientation for 
the best-fit line through the remaining, clustered data are determined. Again, linear regression is 
preferably used to determine the location and orientation of the best-fit line. The negative slope 
(i.e., -1 x m) of the best-fit line provides a pulmonary capillary blood flow measurement, which 
may then be used to determine cardiac output. A correlation coefficient can then be calculated, as 
previously disclosed herein, to indicate the quality of the data used to determine pulmonary 
capillary blood flow or cardiac output. The correlation coefficient or a quality measure based 
thereon may then be communicated to the user (e.g., a doctor, nurse, or respiratory technician) or 
used to weight the resulting pulmonary capillary blood flow or cardiac output value in an output, 
weighted average value. 

[0115] One or a combination of the methods for modifying data may be performed on 
the measured or calculated data to increase the accuracy with which a best-fit line through the 
data or the pulmonary capillary blood flow or cardiac output of a patient can be determined. 

[0116] As an example of the use of filtering and clustering together, the calculated Vco 2 
data are grouped together as the y-axis data of a two-dimensional line graph and the calculated 
CaC0 2 data points are grouped together as x-axis data points. The data points in at least one of 
the groups are filtered to determine a best-fit line for the data having an optimum correlation 
coefficient. The data are also clustered, either before or after filtering, to further improve the 
correlation coefficient of the best-fit line to the calculated VC0 2 and CaC0 2 data. The remaining 
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data is then used to determine (e.g., by linear regression) a best-fit line therefor, as well as a 
correlation coefficient for the best fit line. The slope of the best fit line is then calculated and used 
to determine pulmonary capillary blood flow or cardiac output. The correlation coefficient may 
also be used to indicate the reliability of the pulmonary capillary blood flow or cardiac output 
determination or to impart a specific weight to the pulmonary capillary blood flow or cardiac 
output determination in a weighted average thereof 

[0117] Once the location and orientation of an accurate best-fit line for the data has 
been determined, as disclosed previously herein, the pulmonary capillary blood flow of the patient 
can be calculated as the negative of the slope of the best fit line. 

[0118] In addition, the best-fit line can be used to estimate the CvC0 2 of the patient. 
When VC0 2 eventually ceases during total rebreathing, the partial pressure of carbon dioxide 
measured at the mouth (pC0 2 ) of a patient may represent the CvC0 2 of the patient. When partial 
rebreathing techniques are used, the Vco 2 of the patient is reduced to levels lower than baseline, 
but is not reduced to zero. By determining the best-fit line through data obtained by use of partial 
rebreathing techniques, the best-fit line can be extended to a point where VC0 2 would be equal to 
zero or effectively zero and, thereby, used to determine the carbon dioxide content, or mixed 
venous carbon dioxide content (CvC0 2 ), of the patient's blood at that point. Equation (39), 
which is the equation for the best-fit line, can be rearranged in terms of carbon dioxide elimination 
as follows: 



When carbon dioxide elimination ceases, Vco 2 is equal to zero and equation (50) becomes: 



[0119] Accordingly, the present invention also includes methods for substantially 
noninvasively determining CvC0 2 when partial rebreathing techniques are employed. 

[0120] If the above-described processes are used to determine the pulmonary capillary 
blood flow of a patient, the intrapulmonary shunt flow of the patient or intrapulmonary shunt 
fraction of the cardiac output of the patient may also be determined, as known in the art. The 



Vco 2 = m x CaC0 2 + b. 



(50) 



0 = m x CvC0 2 +b, 
which can be rearranged as follows: 
CvC0 2 = -b/m. 



(51) 



(52) 
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cardiac output of the patient may then be determined from the pulmonary capillary blood flow and 
intrapulmonary shunt flow of the patient, as known in the art. 

[0121] The relatively short phases of dtfj^rential Fick techniques incorporating teachings 
of the present invention, as well as the lack aPa recovery or stabilization period, facilitate the 
calculation and, thus, reporting of noninvasive pulmonary capillary blood flow or cardiac output 
measurements with increased frequejk^y over that possible with previously known differential Fick 
techniques. For example, when conventional partial rebreathing techniques are employed, 
pulmonary capillary blood flow and cardiac output values can only be updated as frequently as the 
cycle time for these methdas, which is typically three minutes or longer. In contrast, when the 
differential Fick meti><xl of the present invention is embodied as a partial rebreathing process with 
rebreathing and n/n-rebreathing phases that last about thirty seconds, the pulmonary capillary 
blood flow and/ardiac output of a patient can be updated following the completion of each 
phase, or abpfut every thirty seconds. 

[0122] While specific rebreathing processes are disclosed herein, the methods disclosed 
herein may also be used with other rebreathing processes, as well as with other differential Fick 
techniques for noninvasively measuring the pulmonary capillary blood flow or cardiac output of a 
patient. 

[0123] The method of the present invention facilitates the measurement of pulmonary 
capillary blood flow or cardiac output on a more frequent basis. As a result, when methods of the 
present invention are employed, the pulmonary capillary blood flow or cardiac output of a patient 
may be better and more accurately tracked. 

[0124] Although the foregoing description contains many specifics, these should not be 
construed as limiting the scope of the present invention, but merely as providing illustrations of 
some exemplary embodiments. Similarly, other embodiments of the invention may be devised 
which do not depart from the spirit or scope of the present invention. Features from different 
embodiments may be employed in combination. The scope of the invention is, therefore, indicated 
and limited only by the appended claims and their legal equivalents, rather than by the foregoing 
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description. All additions, deletions, and modifications to the invention, as disclosed herein, 
which fall within the meaning and scope of the claims are to be embraced thereby. 
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